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We perform an analytical investigation in the framework of generalized K matrix theory of the 
scattering problem in tight isotropic and harmonic waveguides allowing for several open scattering 
channels. The scattering behavior is explored for identical bosons and fermions, as well as for 
distinguishable particles, the main aspect being the confinement-induced resonances (CIR) which 
are attributed to different partial waves. In particular we present the unitarity bounds which emerge 
when considering a quasi one dimensional system. Unitarity bounds are also given for the transition 
coefficients, which show the limitations for efficient transversal (de-)excitations by means of CIRs. 

We analyze the CIR for d-waves and find the intriguing phenomenon of a strong transmission 
suppression in the presence of more than one open channel, which represents an interesting regime 
to be applied in the corresponding many-particle systems. The corresponding channel threshold 
singularities are studied and it is shown that these are solely determined by the symmetry class of 
the partial wave. 

PACS numbers: 34.10.-l-x, 03.75.-b, 34.50.-s 


I. INTRODUCTION 

Trapping ultracold atomic vapors in tight waveguides 
and thus effectively reducing their dimensionality has be¬ 
come a key concept in the contemporary study of ultra¬ 
cold atomic few- and many-body ^sterns, as exotic low¬ 
dimensional quantum phases [IH 3 such as the Tonks- 
Girardeau gas in one dimension (ID) or the Berezinsky- 
Kosterlitz-Thouless transition in two dimensions (2D) are 
available. Besides these intriguing phenomena, the re¬ 
duction of dimensionality also allows for a novel mecha¬ 
nism to control the scattering physics of two-body inter¬ 
actions. This was first pointed out by Olshanii Q who 
studied the influence of a tight cylindrical confinement 
on scattering events. In particular, a resonance appears 
when the s-wave scattering length becomes comparable 
with the length scale associated with the confining waveg¬ 
uide potential. This so called confinement-induced reso¬ 
nance (GIR) as a result can be controlled by tuning the 
trap frequency, which recently also led to the first exper¬ 
imental realization of a super Tonks-Girardeau gas Q. 
A similar prediction of a CIR for spin-polarized fermions 
0 was in the following also experimentally confirmed 0“ 
l9|. Except these experiments discovering CIRs in ID and 
2D, a recent experiment 0 was carried out in mixed- 
dimensions. 

However, the substantial theoretical effort prevail the 
experimental observations, while at the same time build¬ 
ing up a comprehensive understanding of the princi- 
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pies of CIR and suggesting a variety of systems where 
CIRs can emerge. These efforts include works on dif¬ 
ferent waveguide geometries [iM3 , resonant molecule 
formation [l8| , transparency induced by the confinement 
0123, CIRs in mixed dimensions and multiple open 
channels [2ll - [^ . the coupling of various partial waves 
due to the confinement [24, or atomic scattering with 
anisotropic interactions [ 2 ^, [l^]. Alongside with these 
studies we also want to mention recent investigations on 
ultracold quantum gases on atom chips [13-0 where ex¬ 
cited transversal modes of the confinement are utilized in 
order to prepare entangled atom clouds. 

The focus of the present study is the two-body multi- 
mode scattering behavior of atoms in the presence of an 
axially symmetric and harmonic waveguide. In addition 
to the treatment of bosonic and spin-polarized fermionic 
collisions, we also provide a theoretical description of 
the collisional properties of distinguishable particles, for 
which even and odd partial waves contribute. Similar 
to preceding studies 0, 0, the particles are allowed 
to interact with higher partial waves. The constraint 
of colliding with energies below the excitation energy of 
the first excited transversal state is however lifted in this 
work, thus allowing inelastic collisions, where particles 
can be scattered into different channels. The latter are 
asymptotically defined by the transversal trap modes. 
Our approach is based on the fully analytical and non- 
perturbative description in terms of K matrices, whose 
usefulness has already been demonstrated in a series of 
previous works 00, [00. Due to the seminal results 
of Bo Gao 0-0 who investigated the free-space scat¬ 
tering properties of neutral alkaline atoms possessing a 
van der Waals tail, analytical formulas for the (general¬ 
ized) scattering length were derived, which are used in 
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the present setup to adequately describe the scattering 
event on the interatomic scale. 


The K matrix approach provides a generalization of 
the works of Granger et al and Kim et al incor¬ 
porating however all the higher partial waves and contri¬ 
butions from all the closed channels. Furthermore, go¬ 
ing beyond the previous studies we derive the connection 
of the physical K-matrix with all the relevant scatter¬ 
ing observables obtaining thus the full multi-component 
scattering wave function. Using this formalism we study 
the universal properties of ultracold collisions in the i- 
wave single partial wave approximation (GSPWA). Here, 
we find the existence of energies above the correspond¬ 
ing channel threshold at which the collision effectively 
behaves as in free-space. This energy is found to be in¬ 
dependent of the number of open channels for s- and 
p-wave interactions. Furthermore, we investigate the uni- 
tarity bounds for quasi-lD collisions. This quantitatively 
explains the transmission at a CIR in higher transver¬ 
sal modes and explains the confinement induced uni- 
tarity bound, from which we also derive the unitarity 
bounds for inelastic collisions in the waveguide. These 
bounds may be useful when investigating the possibilities 
and limitations of populating higher transversal modes 
by means of a CIR, as interest in coherent excitations 
in waveguides increases Next, we investigate 

the intriguing possibility of a blockade in the first ex¬ 
cited transversal mode, found from the quasi ID unitar¬ 
ity bound in the d-SPWA, with an adequate interatomic 
potential and the coupling of partial waves taken into 
account, showing that an almost totally blockaded trans¬ 
mission channel may exist even if there are other possi¬ 
ble scattering channels available. In addition we discuss 
the scattering of distinguishable particles in waveguides 
and show the qualitative difference of partial wave cou¬ 
pling for distinguishable particles by introducing reso¬ 
nance and transparency coefficients [^. We discuss the 
occurring threshold singularities for collisions of indistin¬ 
guishable particles and show the qualitatively different 
behavior for bosons and fermions. 


In detail, this paper is organized as follows. Section HIl 
gives a brief review of the considered waveguide Hamil¬ 
tonian as well as the applied techniques, namely the K 
matrix formalism while employing the local frame trans¬ 
formation. Thereafter, Sec. imi introduces and discusses 
the relevant scattering observables for ID multichannel 
collisions as well as provides their connection to the phys¬ 
ical K matrix for the case of ID geometries. Section HVl 
is devoted to the analysis of our results while our sum¬ 
mary and conclusions are given in Sec. [V] The Appendix 
provides among others some technical concepts used to 
derive the physical K matrix. 


II. WAVEGUIDE HAMILTONIAN AND 
A-MATRIX APPROACH 

In the following we study the collisional behavior of 
two particles within a harmonic waveguide. Hereby indis¬ 
tinguishable or distinguishable particles are considered. 
The harmonic nature of the confining potential does not 
couple the center of mass (CM) coordinates with the rela¬ 
tive ones and permits us to treat their motions separately. 
The Hamiltonian accounting for the CM motion simply 
describes a CM excitation in a harmonic potential. This 
solution is well known and thus of no further interest. 
The non-trivial part is the relative motion Hamiltonian, 
which reads 

H = + + VLjir), ( 1 ) 

where r = is the interparticle distance, with 2 

and p describing the longitudinal and transversal degrees 
of freedom, respectively, p denotes the reduced mass of 
the colliding pair and u)±_ is the confinement frequency. 
Accordingly, the harmonic oscillator length scale is given 
by a_L = y/h/pujj_. The term ULj(r) = ^ ^ is the 

Lennard-Jones 6-10 potential indicating the short-range, 
two-body interatomic interactions. Cq is the dispersion 
coefficient and it defines the van der Waals length scale 
via the relation Pq = . We regard Cio as 

a parameter in order to tune the corresponding scatter¬ 
ing lengths induced by the short-range potential. Among 
others, the particular choice VLjir) is motivated by the 
fact that we attribute to the two-body physics a realistic 
character avoiding the use of zero-range approximations. 
In addition, this particular type of interactions are an¬ 
alytically solvable by means of the generalized effective 
range theory [^ . However, any other interatomic poten¬ 
tial is also suited as long as the length scale /3„ associated 
with this potential is small compared to the oscillator 
length, i.e. Pn ^ a,±- 

As in previous works [1, [2^ [l^, on the Hamilto¬ 
nian given in Eq. o, the separation of length scales is 
assumed, i.e. Pq ^ a_L. In short, this implies that the 
Hamiltonian has three distinct regions where (a) different 
potential contributions dominate and (b) different sym¬ 
metries are obeyed by the corresponding Hamiltonian. 

(i) Starting in the inner region, where r ^ Pq holds, 
the interatomic potential dominates. The two particles 
thus experience a free-space collision with total energy 
E = 

(ii) The effect of this collision on the wave function is 
best observed from the intermediate region (/Ig <C r <g; 
a_L), where both potential contributions are negligible. 
Hence, we can monitor the outcome of the collisional 
event in region (i) by a well defined phase shift for 
each partial wave. Due to the invariance under rotations 
SO{3) we can arrange the full scattering information in 
a diagonal, energy dependent, K matrix 

(iii) As the asymptotic region r ^ a_L is concerned, 
only the transverse confining potential contributes. The 
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wave function is thus a direct product of a sine or 
cosine function in the z-direction and a 2D harmonic 
oscillator (HO) mode for the p-direction, i.e. I'i/') = 
c„ where denotes the channel momen¬ 

tum in z-direction, n denotes the transversal oscillator 
mode and m the magnetic quantum number. The total 
energy E distributes over these two degrees of freedom 

according to the relation E = huj±{2n -I- |m| -I- 1) -I- 

We note however, that the azimuthal quantum num¬ 
ber m asymptotically associated to the SO{2) rotations, 
is a good quantum number in both regimes and there¬ 
fore is a fixed quantity throughout our analysis which 
is set to m = 0 for what follows and is therefore also 
omitted in the labeling of the states. Hereafter, we drop 
the typical assumption in most of the existing literature 
that the total collision energy has to be sufficiently small 
such that only the energetically lowest transversal state 
can be populated, i.e. |' 0 ) ^ |< 7 o; 0 ) and thus we are 
going beyond previous studies by allowing inelastic col¬ 
lisions involving several transverse modes. Already from 
this expansion which is only invariant under Tz 0 SO(2), 
where Tz denotes translation along the z-direction and 
50(2) rotations around that axis, we conclude that the 
mapping from region (ii) to (iii) cannot be accomplished 
by a unitary transformation. Thus, in order to transfer 
the scattering information between these two regions of 
different symmetry, an appropriate way is given by the 
local frame transformation Uir However, as al¬ 

ready discussed previously [a, the application of this 
technique comes at the price of rendering the closed chan¬ 
nel (fku±(2n + 1) > E) components of the wavefunction 
unphysical. This drawback is due to the boundary con¬ 
ditions of our scattering approach, namely the standing 
wave approach, which after the analytical continuation 
to the closed channels turns the oscillating solutions into 
exponential diverging ones. To overcome this unphysical 
situation a standard channel closing procedure, familiar 
from multichannel quantum defect theory has to be ap¬ 
plied [S^, leading to the physical K matrix, given by 


j^lD,phys _ 


ID 

oo 




( 2 ) 


where in turn, refers to the corresponding ID K 

matrix [I, [13, Si, Si- In addi¬ 
tion, denote the open-open and closed-closed 

channel transitions, respectively. Accordingly, and 
K),^ denote the open-closed channel transitions and vice 
versa. From Eq. the resonant processes are given 
as poles of the physical matrix. Therefore, the 

roots of det(l — iK).^) correspond to the positions of 
the closed channel bound states lying in the continuum 
of the open channels. The resonant structure thus ful¬ 
fills a Fano-Feshbach scenario . Performing the above 
calculations (cf. Appendix [^ for a given 3D K matrix 
_ diag(tantan(5^/) yields a physical K matrix. 


given by 




det(l - *K3^il) ^ 

- - iiw (F^,, + F^,^))}, 

(3) 


where Ag{E) = ta.n6i{E) contains the energy depen¬ 
dent £-th phase shift and the matrices F°g, are given by 
{F°g,)nn' = UinUtn', with 0 < n, n' < Uq - 1 and Uo 
denoting the number of open channels. At this point 
we also introduce the generalized, energy dependent, 
scaled scattering length defined for all partial waves i, 
by ae{E)'^^+'^ = (a^(A)/a_L)^^+^ = - ■ As the 

local frame transformation also depends on the en¬ 
ergy, this also holds for the matrices F°^, = F°^,{E). 
Here, we also note that one should carefully distinguish 
the number of open channels rio from the actual quan¬ 
tum numbers of the corresponding transverse modes n, 
e.g. consider the single mode regime where only the pop¬ 
ulation of the lowest transversal mode |( 7 o; 0 ) is allowed, 
we have Uq = 1 . 

The energy dependent K matrix given in 

Eq. ([3]) appropriately describes the scattering in a tight 
harmonic waveguide with several open transverse modes. 
The coupling of two arbitrary partial waves i and i' be¬ 
longing to the same symmetry class, i.e. £ — £' = 0 mod 2 , 
is properly taken into account. Before proceeding let us 
recall that the effect of the closed channels on the scat¬ 
tering phase shift can conveniently be expressed as 


e+i' 

iltt'(e) = ^ 

p—0 




t^2(e + i)"i 


At 1 „ 


cTn = (_l)^^(2£+l)(2f+ 1) 

e+i' 

X ^ r{£,i\u,p) 

T{t,e = iP-^2''-^{2v + 1) 


X 



(£ i' ie\ 

1^0 0 oj 




(4) 

(5) 

( 6 ) 


where C//(s,a) denotes the Hurwitz zeta function and 
r(-) are some combinatorial constants containing the 
Wigner 3 7 -symbols. Equation (|3]) is discussed in more 
detail in [20j. The dimensionless, channel-normalized en¬ 
ergy e is defined by the relation E = 2fujj±{e + 5 ), which 
is chosen such that n<e<n-|-lis between the thresh¬ 
old of th n-th and the (n + l)-th channel. We note, that 
Uo = [ej -f 1 is the number of open channels, where [ej 
denotes the largest integer smaller than e. The represen¬ 
tation of Eq. o, which differs from the one introduced in 
[ 2 ^, is in particular useful when considering the thresh¬ 
old singularities below. One further remark is in order, 
which refers to the second argument of the Hurwitz zeta 
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function. This argument is given by Ec = Uq — e = 1 — Ae, 
where Ae = e — [ej denotes the fraction of the total colli¬ 
sion energy above the threshold of the last open channel 
rendering (^h periodic with a saw-tooth like behavior. 

Throughout the following analysis the energy ranges 
up to e = 4, passing several channel thresholds of the 
transverse conhnement. We note that all these energies 
lie close to the threshold of the interatomic potential and 
thus the analytic solutions [33| for our interatomic po¬ 
tential are applicable. 


III. SCATTERING OBSERVABLES 


In the asymptotics of the scattering process the 
transversal and longitudinal degrees of freedom are de¬ 
coupled and the quantum number of the transversal 2D 
HO modes can be used to define the asymptotic scatter¬ 
ing channel. 

rio 

Mr) = + E (7) 

n'—O 

describing an incoming wave in channel n which is then 
(in)-elastically scattered into all open channels. Here, the 
scattering amplitude in forward (-I-), respectively 
backward (—) direction reads 

fnn' = fun' + Sgn{z)f °^,, (8) 


whereas in turn and f° refer to the respective scatter¬ 
ing amplitudes for even and odd exchange symmetry and, 
sgn( 2 ;) = zl\z\ denotes the sign function. By the conser¬ 
vation of flux, the forward /"*" and backward f~ scat¬ 
tering amplitude contain the same information about a 
scattering event, we concentrate our analysis on f~^, for 
which the transmission and reflection coefficients Tnn' 
and Rnn’, respectively, from channel n to n' take the 
following form 


= \Snn' + fL' + fL'\^ (9) 

_ \ £p |2 fini 

r^nn' \Jnn' J nn'\ 5 

where i and I' refer to even and odd partial waves, re- 

{i 

spectively. The transition probability characteriz¬ 

ing the transversal excitation and de-excitation processes 
from channel n into a specific channel n' are given by the 
sum of the corresponding transmission and reflection co¬ 
efficients 


’’ n/n.' ^nn.' 


( 11 ) 


If the constituents of the scattering event both belong to 
the same symmetry class, i.e. both are either bosons or 
fermions, we obtain a special case of Eq. given by 
'^nn> = 1'^""' + where the scattering amplitudes 

fnn' connected to the physical K matrix, via 


/W = * K 


lD,phys 

oo,i 



lD,phys 

oo,i 


-1 


( 12 ) 


see also [1^. This relation in particular allows for a 
extension of previous studies on the CIRs in harmonic 
waveguides to distinguishable particles, as it is in detail 
discussed in Sec. IIV Cl 


IV. RESULTS AND DISCUSSION 
A. Universal Properties 

We start our discussion of the universal properties for 
the (in)-elastic scattering in waveguides by considering 
the Uwave single partial wave approximation (USPWA). 
Since the main focus is the study of CIRs which typ¬ 
ically occur in the vicinity of a free space resonance, 
the £-SPWA can safely be assumed to be accurate in 
the description of the scattering process. Recall that 
for the particular choice of a Lennard-Jones type 6-10 
interatomic potential the background scattering length 
from higher lying partial waves {i > 1) is negligible. The 
advantage of the ASPWA is given by the fact, that (i) 
analytical results can be derived straightforwardly, see 
below and (ii) for certain energy regimes it serves as a 
good approximation for the coupled Awave CIR [ISl • 
Employing thus the physical K matrix in the £-SPWA, 
where £ is not restricted to belong to a specific symmetry 
class 


K lD,phys 
00,1 


ai 


, with 


1 — lAtiiu ’ 


(13) 

(14) 


we calculate according to Eqs. ® 
transmission- and reflection coefficients 


^nn' — ^u 


a} {2Ti[¥°u]{ni)un'5n 


l-«?Tr[F?,] 


and OT the 

- jFMun') 


R 


W 

nn' 


l-«£"Tr[Fa 2 ’ 


(15) 

(16) 


from which we derive in particular the total transmission 

coefficient Tn'^ = J2n'<n '^in'^ when incident in channel 
n which takes the appealing form 


tP = 1 - 


'T£\X^^r>,physyr^lD,phys^ 


1 -f 


(17) 


This quantity is in the focus of the following analysis since 
it encapsulate all the relevant scattering information for 
both all the open and closed channels. 


1. General aspects of eollisions in harmonic waveguides 

In this subsection we focus on the general behavior 
of the transmission coefficients for partial waves £ = 0 
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and £= 1 at total collision energies beyond the single 
mode regime. In particular we note, that the separation 
of length scale induces an energy scale separation. This 
means that even several quanta of the transversal excita¬ 
tion imply that the corresponding energy dependence of 
the scaled s-wave scattering length is negligible. Due to 
the increasingly narrow width of higher partial wave res¬ 
onances, this simplification is not legitimate for £ > 1. In 
addition we should recall that all the scattering lengths 
used in the analysis below are analytically obtained via 
a Lennard-Jones 6-10 potential. 



(ai)^ 


FIG. 1. (Color online) (a) Transmission coefficients 
versus the scaled s-wave scattering length no, 
for the hrst, second and third open channel (solid, dashed and 
dotted lines), respectively, at energy Ae = 0.95. Higher lying 
curves for a particular number of open channels correspond 
to a lower entrance channel. Panel (b) shows 
at Ae = 0.05 for the case £ = 1 . 

First we show a typical case of the transmission coef- 
ficients T„ for £ = 0,1 versus the corresponding scaled 
scattering length in panels (a) and (b) of FiglU 

respectively. In particular we observe the asymmetric 
line shape of the transmission spectra for £ = 0. The 
scaled energy above the corresponding channel thresh¬ 
old Ae = 0.95 is chosen such that the CIRs, identified 
as the minima of the respective transmission coefficients 
are best pronounced, i.e. we aim at a large difference 


between the transmission values taken for large |do| and 
the specific value of uq leading to a CIR. This relative 
difference maximizes especially at energies below every 
threshold which can be readily seen in panel (a) of Fig. 
[5l where the transmission coefficient values for |do| —?► oo 
(T^\ - black lines) and for do at a CIR (T^jj^ j - red lines) 
are displayed. Furthermore, we observe that a transmis¬ 
sion blockade is present only for the single mode regime 
(see Fig[T] (a) red solid line). However, in the case of 
several open channels the transmission blockade is lifted 
giving in turn rise to finite values of the transmission co¬ 
efficient. In addition transparency, i.e. = I occurs in 
the absence of interactions = 0) between the col¬ 

liding particles. Analogously, panel (b) of Fig. [T] depicts 
the results for £ = 1, as a function of the scaled p-wave 
scattering volume (di)^. Unlike the s-wave case, here 
the asymmetry of the line shape is barely visible if the 
collision energy is raised above the first excited channel 
threshold. 



FIG. 2. (Golor online) (a) Transmission coefficient versus 
the channel scaled energy e for £ = 0 for various scaled s- 
wave scattering lengths. The GIR obeys = 0 in the 
ground channel and for higher channels the resonances follow 
the confinement induced unitarity bound q (thick, dot- 

dashed line), see Eq. (1221) . (b) versus e for the same 
scattering lengths. Again, the resonances are bounded by 

7^(0) 

More insight into the behaviour of the transmission 
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coefficients can be obtained from Fig. [2l where we 
present and versus e in panels (a) and (b), re¬ 
spectively. These results correspond to different values 
of the Cio parameter. As it is discussed in more detail 
below (see Sec. IV A subsection 3 and 4), a CIR process 
belongs to the transmission minima which occur in the 
interval of energies among the open channels. Specifi¬ 
cally, in Fig. [2] (a) we observe that for e < 1, ie between 
the thresholds of the ground and first excited transverse 
state, the transmission minima for varying values of the 
scattering lengths leads to a transmission blockade. For 
1 < e < 2 and higher energies the minimal value of 
the transmission is nonzero and increases with increas¬ 
ing channel the energy belongs to. The transmission ex¬ 
hibits a repeating pattern shifted to higher values. The 
locations of the minima with varying scattering length 
form thus a “topos“ which in Fig. [2] is denoted by the 
black dot-dashed line: This is the confinement-induced 
(Cl) unitarity bound and represents a universal feature. 
We should note however that the specific functional form 
of the Cl unitarity bound depends on the particular £- 
wave character of the collisions. Furthermore, in Fig. 
Ha) we observe that a CIR occurs for positive scattering 
length at Ae < 0.69, while CIRs of negative values of oq 
emerge for Ae > 0.69. The particular value of oq = 0.3 
does not permit a CIR at all and is hence monotonically 
increasing, while departing from e = 0 with a finite, non¬ 
vanishing slope. 

Finally we note that the transmission coefficient goes 
up to unity at the channel thresholds regardless of the 
value of the scattering length. Considering the case ^ = 1, 
presented in Fig. H we use the Cio parameter to label 
the different curves. This is due to the narrow width of 
the p-wave free space resonances which makes an energy 
independent treatment of the scaled p-wave scattering 
length impossible. For each curve the Cio parameter is 
adjusted such that a free-space resonance occurs for every 
energetic interval i — 1 < e < i for 0 < i < 4. Again, at 
a CIR the corresponding transmission coefficient touches 
the p-wave Cl unitarity bound and thus we encounter a 
suppression of the transmission with a complete block¬ 
ade T = 0 for the case e < 1. Away from a free-space 
resonance the scattering length quickly decreases to its 
small background value leading to the large value of the 
transmission coefficient away from a resonance. Further¬ 
more we observe that in the case of p-wave interactions 
the value of taken at the channel thresholds indeed 
strongly depends on the scaled p-wave scattering length, 
which drastically differs from the case of £ = 0 thus ren¬ 
dering the threshold behavior for £ = 0 universal with 
respect to the s-wave scattering length. 

2. The decoupling energies 

In Ref. [1^ it was shown that the coupling of the par¬ 
tial waves to the closed channels, given by the il«(e)’s in 
Eq.dD) exhibits roots for all considered partial waves £. 



FIG. 3. (Color online) (a) versus the scaled energy e for 
four open channels. Different curves correspond to different 
values of the short range parameter Cio. It is clearly observed 
that the CIR saturates at the confinement induced unitarity 
bound depicted as the dot-dashed curve. As for general i a 
CIR associated with a transmission blockade is only happen¬ 
ing for e < 1. In (b) the transmission coefficient Ti^\ when 
incident in the first excited channel is shown. 

These roots -b Ae|^, in the following called 

decoupling energies, depend on the partial wave £, and, in 
general, also on the number of open channels n. The de¬ 
coupling energies determine the particular energy where 
the bound state of the closed channels decouples from all 
the open channels. Therefore, in this case, the physical 
K matrix can be written in the following form 

= (K3^)« F^,(LeJ+e*) 

= U^K3^U|l.j+.. (18) 

where [ej is the threshold energy for the largest open 
channel. Equation (TTSl) shows the expected result, that 
the 3D scattering information which emerges close to the 
origin is transfered to the asymptotic regime without be¬ 
ing affected by the closed channels of the trapping po¬ 
tential. In other words the colliding pair experiences 
effectively a free space collision in the presence of the 
waveguide geometry. 

In particular by using the expression for ilu{e) in Eq. 
(H one can show, that for the cases £ = 0 and I, the de¬ 
coupling energies do not depend on the number of open 
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FIG. 4. (Color online) (a) Transmission coefficients at the 
decoupling energies [eoj + ^^0 versus the scaled s-wave scat¬ 
tering length ao. Different curves depict the transmission co¬ 
efficients for up to three open channels. The symmetric line 
shape centered around the non-interacting case ao = 0 in¬ 
dicates an effective decoupling from the closed channels of 
the waveguide. For large values of |ao|, the transmission ap¬ 
proaches the values given by Eq. dsni). The case of ^ = 1 is 
shown in (b) for energies [eij -I- Ae*. 


channels, i.e. = Ae|, which means that the closed 

channels decouple from the open channels at the same 
energy Ae| above the channel threshold, as it is illus¬ 
trated in Fig. m Fig. SKa) shows the bosonic case, with 
i = 0 and Acq « 0.69, where the transmission coefficients 
jj(o) ^ j,(o) ^ j,(o) .^gj-g^g scaled s-wave scattering length 
do are shown for one, two and three open channels, re¬ 
spectively. The symmetric line shape centered around 
the non-interacting case do = 0 is clearly seen. Trans¬ 
missions for £ = 1 and e* ~ 0.31, are shown in panel 
(b) as a function of the scaled p-wave scattering volume 
(di)^. In both panels the respective scaled scattering 
length are considered as external parameters which is in 
the £-SPWA equivalent to a change of the Cio parameter. 


3. The CIR limit 

As already emphasized above a CIR occurs in the vicin¬ 
ity of a free-space resonance, i.e. in parameter regions 
|d|| 1. For this free-space unitarity regime the scat¬ 

tering amplitude matrix in the CSPWA is simply given 
by 


T70 

M{i) ^_ 

iI« + Tr[F,°]’ 


(19) 


from which we readily derive the corresponding transmis¬ 
sion coefficients 


'tW 

00,71 


- Tr[F° ](Tr[F° ] - (F° )„„) 
K - Tr[F?,]2 


( 20 ) 


This equation gives the values in the wings of large scat¬ 
tering length of the transmission coefficients of Figs. [T] 
and m Furthermore, its general energy dependence is 
show in panels (a) and (b) of Fig. [5] for the cases of 
£ = 0 and 1, respectively. There the transmission coeffi¬ 
cients (ie. T^]n) for large values of the scattering length 
are compared to the quasi ID unitarity bound, i.e. the 
bound to the transmission coefficient at a CIR ie. Tcir 

To understand the behavior of the transmission coeffi¬ 
cient at a CIR, recall that the denominator of Eq. m 
represents det(l — iK),^), which implies, that if this ex¬ 
pression vanishes for a particular value of A^, or the cor¬ 
responding scattering length, respectively, a CIR occurs. 
We hence have a sufficient criterion for the occurrence of 
a CIR given by the condition —?► 00 , with ae from Eq. 

(ED- A transmission value of T = 0 can only be achieved 
in the s-SPWA when there is a single open channel, since 
then Tr[F°^] = F°g holds. Erom a physical point of view 
this behavior is expected, since a transmission blockade 
in a specific channel will not prevent the scattering into 
other channels, which for example can be seen in Eig. 
[3 where the transmission coefficients and are 
shown, respectively. Both curves exhibit minima at the 
CIR. Also in Eig. [7] we plot the transition coefficient 
(dotted lines) which is enhanced at a CIR clearly 
demonstrating that for two open channels at CIR the in¬ 
elastic process are enhanced prohibiting in this manner 
the transmission coefficient to be zero when more than 
one channel is involved. 

The absence of a transmission blockade at a CIR in 
the case of more than one open channel, can be derived 
quantitatively from the formal limit —>■ 00 taken in 
Eq. (fTTll . This yields 


fW __ 

•/CIR ~ 


F o 


Tr[F^ 


( 21 ) 


which generalizes the well known single mode result 
/ciR ~ It is remarkable that for the transver¬ 

sal ground state a non-trivial energy dependence is not 
present. We note that Eq. m can, by virtue of [^ . 








































FIG. 5. (Color online) (a) Transmission coefficient ^ at 
a CIR, i.e. the confinement induced unitarity bound, versus 
the scaled energy e as well as the transmission coefficient 
taken for large scattering length |a^|, i.e. for scattering at the 
free-space unitarity bound, (b) Corresponding result for £ = 
1. Note that the two curves rapidly converge, implying a far 
less pronounced CIR in terms of a transmission suppression 
than in the case of £ = 0. 

also be obtained by formally equating Eq. m at a van¬ 
ishing closed channel coupling, i.e. fci^ — /oo^|um=o- 
Writing down the transmission coefficients for the scat¬ 
tering amplitude at a CIR, or, equivalently considering 
^cm,n = yields 

rpW _ 1 _ {^u)nn 

CIR." Tr[F^J ■ ^ ' 


sion. The difference between the transmission for infinite 
scattering length and at a CIR is largest close to the 
channel thresholds for ^ = 0. Remarkably, this difference 
^CiR i ~ for ^ = 1 happens to vanish as soon as the 
channel threshold to the first excited mode is exceeded. 
This also explains the rather symmetric line shapes ob¬ 
served in panel (b) of Fig. [T]for Uo > 2. 



e 


FIG. 6. (Color online) Confinement induced unitary bounds 
^cmnn' foi' energies e < 3. Panels (a)-(c) depict the cases 

(i) 

£ = 0,1 and 2, respectively. In particular, coincides 

for £ = 1 and 0 for the first excited channel, as can be seen 
in panels (a) and (b). The d-wave bound (c) exhibits a non¬ 
monotonic behavior as well as a total suppression of transi¬ 
tions at e = 1.75. 

A similar analysis can be performed in the ASPWA 
for the transition coefficient In this case, the limit 

ttf —>■ oo leads to 


This expression in particular contains the previous state¬ 
ment that a CIR for a single open channel has = 0. 

In the ^-SPWA the smallest value the transmission co¬ 
efficient can take as a function of the energy is determined 
by Eq. ((22ll . which serves as a lower bound for the trans¬ 
mission coefficient. This lower bound is compared to the 
value of the transmission coefficient for infinite scattering 
length in Fig. [5] The channel of incidence varies within 
the first few transversal modes, i.e. 0 < i < 4. As it can 
be easily deduced from Eq. (1^ all quasi ID unitarity 
bounds tend to unity irrespective of the partial wave £ 
under consideration. This makes the CIR being less pro¬ 
nounced for increasing energies (excited channels) with 
respect to the corresponding suppression of the transmis¬ 


Tjr(^) _ o ££')nn' 

CIR.nn' - ^ Tr[F°J2 ’ 


(23) 


describing the unitarity bounds for the transition coef¬ 
ficients which occur between different channels during a 
CIR. Results for this observable are shown in Fig. [51 
where all possible transition coefficients are given for en¬ 
ergies e < 3. In particular, panel (a) and (b) depict 
the cases for £ = 0 and 1, respectively. In the ener¬ 
getic range 1 < e < 2 the transition coefficients exactly 
coincide and exhibit a monotonically increasing behav¬ 
ior. This coincidence abruptly changes when crossing the 
channel threshold to the second open channel especially 
for VEcir 121 which has a sharp drop to zero for £ = 0 
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while for £ = 1 the coefficient continuously decreases but 
remains the dominant transition process throughout the 
channel. Figl6|c), depicts the case for £ = 2, which ex¬ 
hibits a non-monotonic behavior for all considered transi¬ 
tion coefficients. Even though, as discussed in more detail 
below in Sec. IIVB 11 we do not expect this behavior to 
exactly describe the processes for d-waves as additional 
s-wave contributions have to be taken into account, it is 
nevertheless remarkable that for £ = 2 energetic regions 
exist where by means of a CIR no transitions between 
channels can be induced, i.e. at the energy e = 1.75 the 
transition probability between the two available channels 
vanishes. 

The interest in the observable W is given by recent 
experiments on atom chips [27l - [^ where coherent ex¬ 
citations in higher transversal confinement modes are 
engineered. From this viewpoint the presented analy¬ 
sis may contribute to an understanding in how far the 
CIR may be utilized to coherently excite atoms to higher 
modes and which are the most efficient energetic regions 
in which these transitions can be achieved. 


4- Unitarity and CIR 

To further illuminate the unitarity regime let us con¬ 
sider the situation from the viewpoint of traditional scat¬ 
tering theory. Here, by using the appropriate relation for 
one spatial dimension 


S = l + 2f (24) 

between the scattering matrix and amplitude we 
readily derive the quasi ID unitarity relation 

(25) 

where the right hand side denotes the real part of the 
scattering amplitude matrix, whereas the left hand side 
in particular encapsulates the total reflection coefficient 

= J2n<no^nl' = iff^)nn- We thus Con¬ 

clude from the unitarity relation, that the total reflection 
coefficient when incident in a specific channel n is 
fully contained within a single element of the scattering 
amplitude, namely 

RW=-91e(/W) (26) 


Equating now the right hand side with our system specific 
information from Eq. (ini), we find after some algebra 



{F°g)nn 1 

1 -I- 


(27) 


which shows, that the scattering saturates at the uni¬ 
tarity bound ae ^ oo. It is hence legitimate to regard 
^CiRn equivalent quantity Rn'^\af^oo as the con¬ 

finement induced (Cl) unitarity bound, similar to the 
unitarity bound in free space, which scales as k~'^. 


B. Transmission suppression in excited channels 



FIG. 7. (Color online) Transmission coefficient (red solid 
line) for a small negative value of the s-wave scattering length 
where the CIR is expected to occur near the threshold, namely 
at £ = 1.98. At the same position occurs also the CIR in the 
first excited channel see (dashed line). The total blockade 
is absent, since the transition (dotted line) between the 

channels is also resonantly enhanced at the position of the 
CIR. 

Increasing the total collision energy across the thresh¬ 
old of the first excited transversal mode lifts in general 
the blockade in a particular channel as the scattering con¬ 
stituents may escape via inelastic collisions to a different 
asymptotic state. This is for example present in Eig. 
[71 where a pronounced s-wave CIR is shown for a total 
collision energy of e ~ 1.98. In particular, the transmis¬ 
sion coefficients and are shown together with 
the transition amplitude illustrating the resonantly 
enhanced transition between the channels at a CIR. 

Although Fig. |7]depicts the case of s-wave interactions 
it provides the typical scenario when a CIR happens in 
the presence of more than one open channel. However, 
the first exception where there is no resonantly enhanced 
transition between the channels at a CIR, occurs for d- 
waves when the particles are incident in the first excited 
channel as we will demonstrate in Sec. IIV B 21 But before 
doing so, we will provide a general brief discussion of the 
influence of higher partial waves. 


1. The impact of higher partial waves 

The analysis of s- and p-waves scattering carried out in 
the £-SPWA was justified because of the negligible phase 
shifts of allowed higher partial waves in the presence of 
a free-space resonance of £ = 0 and 1, respectively. This 
particular simplification does of course depend on the 
considered interatomic potential, and in case of Lennard- 
Jones 6-10 potential this certainly holds. Nevertheless, 
when considering the CIRs associated to a higher partial 
wave, e.g. d-wave for the present discussion, we have to 
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FIG. 8. (a) versus e for coupled partial waves I = Q and 

2 constituted by different regions of near resonant behaviour 
belonging to different Cio values. In the middle of each verti¬ 
cal stripe a free-space resonance translates to a corresponding 
CIR. The thick dashed line depicts q. Deviations of the 
transmission minima from the unitarity bound arise due to 
the d-SPWA. The inset (b) shows a zoom-in plot around the 
blockade region [T = 0) in the neighbourhood of e « 1.75. 
The leftmost resonance (see arrow) coincides with the uni¬ 
tarity bound in which case the partial waves decouple, i.e. 
il 20 = 0. Panels (c) and (d) show the results for Devia¬ 

tions from the d-SPWA become small around the decoupling 
energy. We also observe that around this energy scattering 
becomes almost transparent. 


take into account the non-vanishing s-wave scattering in 
the case of collisions between indistinguishable bosons. 
This observation led to the idea of partial wave coupling 
due to closed channels of the confinement, as firstly dis¬ 
cussed in [^. The corresponding physical K matrix of 
this problem based on Eq. ([3]) is now expanded in terms 
of r°^, F°^, and F°,^, which serve as a basis for the 

two partial wave K matrices. Therefore, it is not feasible 
to obtain a version of Eq. (HU) valid for two partial waves, 
since an inversion as needed in Eq. m by means of 
a repeated application of the Sherman-Morrison method 
would yield a physical K matrix with several hundreds of 
summands and is thus prohibitive. We therefore have to 
rely on a more qualitative analysis and take the intuition 
from the transparent results obtained in the £-SPWA. 

However, the coupling between the two contributing 
partial waves i and mediated by the closed channels, is 
conveniently described by the corresponding off-diagonal 
element of the closed channel coupling matrix it. 

Similar to the previous case in Sec. IIV A 21 where the 


vanishing diagonal elements it^^ were used to identify the 
decoupling energies for which the scattering process ex¬ 
hibits strong free-space character within the confinement, 
we generalize it here in the presence of a i'-i wave sys¬ 
tem. Our first observation is that there exist also en¬ 
ergies for which the off-diagonal elements of it vanish, 
i.e. = 0 and the confinement-induced coupling 

between the partial waves is absent. For the case of a 
coupled s-d-wave system, this energy in the first excited 
channel is eo ~ 1.65. In particular this also implies that 
the position of the d-wave CIR is solely determined by 
the d-SPWA, i.e. 


RC2o(eD) = RC2(eD), (28) 

where the resonance coefficients RCf(e) and RC«'(e) [2^ 
are described in Appendix [B1 

For the case of s- or p-waves, the coupling to the closed 
channels has a discrete shift symmetry, i.e. 

il«(e -I- IV) = il«(e), with NgN, (29) 

which is equivalent to saying that the asymptotically de¬ 
fined transversal modes have no influence on the coupling 
through the closed channels in the lowest possible partial 
wave for each respective class of exchange symmetry. The 
origin of this symmetry is clearly the equidistant modes of 
the transversal harmonic oscillator in combination with 
the particularly simple nodal structure of the local frame 
transformation for s- and p-wave interaction for the case 
m = 0. This changes from the d-wave on. Since this fea¬ 
ture does not depend on the coupling of partial waves, it 
can already be understood in the d-SPWA, where a CIR 
occurs when d^ = RC^(e) is fulfilled. By inserting Eq. ([4]) 
in Eq. jui) one obtains 


df+i 


2e ,_ 

P—0 


p-1 
2 


Cc 




(30) 

which depends on the total collision energy e for partial 
waves other than i = Q and ^ = 1, since = 0. 

This rather technical observation of a channel dependent 
coupling to the closed channels may indeed become rel¬ 
evant when experimentally trying to observe a d-wave 
CIR for collision energies e > 1, since the correspond¬ 
ing scattering length required to be comparable with the 
confinement length scale is reduced by the additional en¬ 
ergy dependent factor in Eq. (1301) . which will likely make 
the d-wave CIR less difficult to be observed at energies 
e > I. Having clarified this we can now focus our dis¬ 
cussion of the transmission suppression in the regime of 
multiple open channels. 


2. Strong suppression in the first excited channel 

Let us start the discussion on the strong transmission 
suppression in the first excited channel by a qualitative 
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analysis of the d-SPWA from which we will gain some 
insight. A sufficient and necessary condition for the oc¬ 
currence of a transmission blockade in the first excited 
channel (1 < e < 2) is of course a vanishing of the el- 
ement for either n = 0 or 1, respectively. In¬ 
deed this happens in the d-SPWA for at an en- 

ergy et = 1.75. At the same energy i has to acquire 
an extremal value equal to unity. These two observations 
imply that a CIR at ei, will lead to a blockade for particles 
which are incident in the lowest transversal mode, while 
for particles incident in the first excited channel, the CIR 
will result in complete transparency, i.e. 1 = 1- 

Taking additionally into account s-wave scattering 
which couples to the d-wave this will influence the above¬ 
observed phenomena. The corresponding results are pre¬ 
sented in panels (a) and (c) of Fig. [8j where the two 
transmission coefficients and are shown ver¬ 

sus e and the notation 2g is used here to avoid confusion 
with Eq. ([2]) but to emphasize that the background con¬ 
tribution from the s-wave is properly taken into account. 
Here, the individual resonances lying within the shaded 
vertical stripes result from changing the Cig parameter. 
Due to the narrow width of the d-wave resonance only a 
relatively small window on the energy axis is relevant to 
the collisional process since the scattering length rapidly 
decays back to a very small value, resulting in a trans¬ 
mission close to unity. From these two panels one readily 
observes that the d-SPWA for the Cl unitarity bound 
^CiRO/i’ is shown by the black dashed curve, is 

a rough approximation to the coupled d-s-wave system. 
However, by probing the Cl unitarity bound for the cou¬ 
pled system we see that this approximation qualitatively 
captures the observed values of the transmission coeffi¬ 
cients at a CIR. 

Furthermore, we observe in panels (b) and (d) of the 
same figure, that around Cf, the d-SPWA becomes more 
reliable, which may be related to the decoupling of the 
partial waves, i.e. ilg2 ~ 0 also for energies in that re¬ 
gion. Nevertheless, the blockade expected in channel one 
from the d-SPWA, as well as the transparency in channel 
two are not observed. This however is totally expected 
since the non-negligible s-wave scattering length in the 
background prevents these extremal values. But still, 
for an extended region of energies 1.6 < e < 1.8, the 
portion transmitted when incident in the ground chan¬ 
nel is less than 1%. Similar, the transmitted part when 
incident in the first excited channel is above 99%. Fur¬ 
thermore, the vanishing of the transmission coefficient in 
the ground channel, i.e. = Tgg"^ + Tqi°^ « 0 in par- 

/q 

ticular also implies that Tq^° rs 0 , which in turn leads 
to VFgi“^ ~ 0. This means that in this region of ener¬ 
gies the scattering processes preserve the channels, i.e. 
elastic collisions dominate since the transition probabil¬ 
ity between the two open channels is negligible. We em¬ 
phasize that this behavior corresponds to a blockade for 
particles in the ground state while excited particles are 


effectively non-interacting regardless the fact that there 
are two open channels. This observation might lead to in¬ 
teresting implications for the corresponding many-body 
system, like a mixture of a Tonks-Girardeau gas and a 
non-interacting gas. 


C. Distinguishable particles 


Let us start by briefly stating some single channel re¬ 
sults, for which a similar analysis as in the case of indis¬ 
tinguishable particles was carried out before. As intro¬ 
duced in Eq. ®, the appropriate scattering observable 
for distinguishable particles (DP) is the transmission co¬ 
efficient where both scattering amplitudes, one even and 
one odd, are present. The fact that the Hamiltonian of 
distinguishable particles still commutes with the parity 
operator permits us to treat collisional events within the 
framework of the K-matrix approach. Therefore, here¬ 
after we will employ the Ansatz from Eq. (jH]) with two 
single partial waves, namely s- and p-wave, yielding 


nn' 


\ ^ / nn' 


nn' 

(31) 


and denote the corresponding K ma¬ 

trices for the s- and p-wave, respectively. We note, that 
likewise also both K matrices with two contributing par¬ 
tial waves could be used, to describe a system where s-, 
p-, d- and /-waves are significant. From the scalar, i.e. 
single channel version of Eq. an effective physical 

K matrix can be constructed which is explicitly given by 


T^lD,phys 

'^oo,eff 


n^lD,phys T^lD,phys 

^oo,s ^oo,p 

-L “T -tiOO,S -^00,P 


(32) 


The construction is accomplished by simply demanding 
that the transmission coefficient takes the usual form of 
T = (1 -I- [25. As in our previous study, this 

effective K matrix is clearly separated in numerator and 
denominator, which gives rise to transparency and res¬ 
onance coefficients TC«/ and RC«', respectively. These 
coefficients are given by 


RC()';^(e) = - - - 




(33) 


Tc!y(«) = ' 


2^,le+\ 


21'+ 1 


ApUfi 


C/2g - ’ 

(34) 


where the superscript (dp) abbreviates “distinguishable 
particles” to separate the notation from the one used in 
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Eqs. (IB4I) and (IB6I) for indistinguishable particles. Here, 
I and a belong to different, namely even and odd, sym¬ 
metries. The explicit form of these coefficients is very 
similar to the ones introduced in (2^ but, unlike them, do 
not contain off-diagonal elements of the matrix il, which 
indicates the absence of couplings between even and odd 
partial waves due to the confinement, becoming evident 
when comparing Eq. HMD from above with Eq. (lB4l) . 
corresponding to indistinguishable particles. If we inter¬ 
pret Ug'oUio as the trace over the open channels, which 
indeed can be shown rigorously in a lengthy calculation 
for a problem involving multiple open channels, we see 
that the term proportional to a( in Eq. (l33l) is this trace 
squared, contrasting Eq. (IB4I1 where a analog squared 
trace over the closed channels is present. The same ar¬ 
gument also holds for the transparency coefficients Eqs. 
(l34)l and (IB6I) . showing that the mechanism relating the 
partial waves is different for indistinguishable and distin¬ 
guishable particles, respectively. Whereas the former are 
coupled through the closed channels, the latter are con¬ 
nected via the frame transformation in the open channels. 
Even though, there is an obvious structural similarity 
between the corresponding resonance and transparency 
coefhcients for distinguishable and indistinguishable par¬ 
ticles the physical mechanism coupling the partial waves 
is very different. Since the parity in ^-direction is a good 
quantum number the local frame transformation is not 
allowed to couple partial waves with A£± 1. The appear¬ 
ing “mixing term” proportional to ai in the coefhcients 
and originates thus from the coherent su¬ 

perposition of the partial s- and p-wave contributions 
assumed in Eq. (IMJ. 

The arguments for the single channel case made above 
are illustrated in Fig. [9l Here, the black solid line shows 
the transmission coefficient versus Vq, the depth 

of a spherical square well potential which was used in 
order to model the two-body interactions of the DP. For 
comparison there are also the corresponding transmission 
coefhcients shown for the pure s- (red dashed) and p-wave 
(blue dotted) interaction, respectively. It is clearly ob¬ 
served that a dual CIR corresponding to complete trans¬ 
parency appears when the transmission coefh¬ 

cients of s and p are equal, i.e. by Eq. (15^ this means 
that the quasi ID wavefunctions induced by s and p wave 
interaction possess phase shifts of equal magnitude but 
differ by a sign yielding thus destructive interference. As 
in [l9j |. this potential is used to mimic the possibility of 
having large values of Og and Up simultaneously, illustrat¬ 
ing the peculiar quasi ID feature of total transparency 
while strongly interacting with two partial waves. How¬ 
ever, by inspecting the resonance and transparency coef¬ 
hcients and we see, that it is not necessary 

to have large values for both scattering lengths, i.e. a 
sufficient condition is one of them being large and the 
other one possessing a small but non-negligible value, as 
it was already discussed in [2^ for indistinguishable par¬ 
ticles. In addition, the inset shows the intersection of the 
s- and p-wave scattering length depicted by the red and 



Vo of a spherical square well. Black solid lines: with s- and 
p-wave interactions; red dashed line: s-wave approximation; 
blue dotted line: p-wave interactions. The dual CIR is clearly 
seen to happen where the individual transmission coefficients 
of s- and p-wave interactions are equal, as well as the shift 
of the p-wave CIR due to the non-negligible s-wave scatter¬ 
ing. The inset shows the scaled s- (red solid) and p-wave 
(blue solid) scattering lengths versus the potential depth. The 
crossings with the DP transparency coefficients (dashed lines) 
determines the position of the dual CIR. 

blue solid lines with the corresponding transparency coef¬ 
ficients TCjg’’^ (red dashed line) and TCq'^'’^ (blue dashed 
line), respectively. 

Let us now address the multichannel scattering proper¬ 
ties, based on Eg. (1511) . We note, that an expanded repre¬ 
sentation for multiple open channels, similar to Eq. CH), 
in the single mode regime, is straightforwardly derived 
and thus not presented here. By considering a realistic 
short range potential as before, i.e. Lennard-Jones 6 — 10 
we have to distinguish two cases. Firstly, the s-wave CIR 
in the presence of p-wave interactions and vice versa. 

An inspection of the corresponding scattering lengths 
in the former case shows that the p-wave scattering length 
in the vicinity of a s-wave free space resonance is negli¬ 
gible, as it is the case for the indistinguishable s — d- 
coupled case. Therefore, also for distinguishable parti¬ 
cles, RC[,f^(e) ft! RC[,^P^(e) := RCo(e), i.e. the s-SPWA 
is a very good approximation. Hence, the transmission 
coefficient versus energy would result in a behavior very 
much like the one shown in Fig. [2] However, we note, 
that this simplification follows from our choice of L J-type 
potential, which, as already said before, is the adequate 
choice to describe the interactions of neutral atoms. Nev¬ 
ertheless, in the vicinity of a free-space p-wave resonance 
the s-wave scattering length possesses a non-negligible 
value and thus its contribution has to be taken into ac¬ 
count. The result on the transmission is shown in Fig. [TUI 
where we consider the cases of incidence in the two lowest 
channels for the case of four open channels. To obtain 
the different curves, the Cio parameter is adjusted, such 
that there is a free-space resonance of p-wave character 
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FIG. 10. (a) Transmission coefficient versus the scaled 

energy e (b) transmission coefficient when incident in 

the first excited channel. Similar to Fig. [3] individual curves 
correspond to different choices of the Gio parameter in order 
to have resonanes in every energetic interval i — 4 < e < i. In 
addition we show in both panels (a) and (b) the corresponding 
Cl unitarity bounds ^ (black dotted curve) for n = 0,1. 


at a particular energy. Focusing on the case of p-wave 
CIRs for distinguishable particles, the black dotted line 
represents Tcjr i in Fig. 10(a), while 2 depicted 
in Fig.10(b), showing that the neglect of a background s- 
wave scattering length leads to increasing deviations for 
more open channels from the unitarity bound. However, 
we generally observe that the transmission is a mixture 
of s- and p-wave scattering. For the s—wave part this is 
best seen by the zero energy value Ti = 0 and the e = 1 
value of T 2 which is also zero. For the p—wave part we 
observe a transmission coefficient at threshold which is 
less than unity and strongly depends on the Cio param¬ 
eter. These are two features we find to be present in 
the p-SPWA. Therefore the contributions from both ex¬ 
change symmetries strongly contribute to the scattering 
physics of distinguishable particles. 


D. Threshold singularities 

In Figs. El [5] and [TUI we observe that at every channel 
threshold the transmission spectra exhibit kinks. This 
constitutes another aspect of inelastic collisions, the so- 


called threshold singularities. More specifically, it is 
known [4l|, |4^ that when the total collision energy leads 
to the opening of a new channel where new states be¬ 
come available this leads to a non analytic behavior of 
the scattering matrix elements. 

In order to firmly address this point it is useful to first 
consider the behavior of the individual elements consti¬ 
tuting the AT matrix when energetically approaching a 
channel threshold from below or above the closed chan¬ 
nel thresholds. Since the notion of transversal channels 
is inherently connected in our framework with the lo¬ 
cal frame transformation @ and its derived quantity 
ilf£', the channel threshold behavior will solely depend 
on these quantities and their properties around a corre¬ 
sponding threshold. We start this analysis with the local 
frame transformation [//„, for which we find the following 
expression for the limit from below 


(-i)"°V2rTT 

e% [4(iV-m)(fV-bl/2)]i/4^'^y Ar + 1/2''’ 

(35) 

which holds for all open channels m < N, where N 
denotes the threshold to the lowermost closed channel. 
From Eq. (1551) we observe, that the elements of the lo¬ 
cal frame transformation acquire a finite value at thresh¬ 
old and are continuous across the threshold. Next we 
inspect the limit to the channel threshold from above 
for the element of the local frame transformation which 
corresponds to the least open channel, i.e. UiN, which 
will then become a closed channel when slightly further 
decreasing the energy e. For this element we find the 
following expressions, which have to be distinguished for 
even and odd partial waves 


Jta C/„ = C+„ + 0{A4)), 


lim UiN = 
t\N Ae\,0V 


^IN — 


lim + 0{Aei 

^2/ 


1 n 


[(Z/2)!]22'+^(iV+i)i/4 


1 (T 


2^ 

1 )! 




with 


(36) 

(37) 


where Eqs. (IM)) and dSll) refer to the case of even and 
odd partial waves, respectively. and denote 

constants depending on the number of open channels as 
well as on the partial wave. Both equations exhibit a 
singular behavior at threshold when approaching from 
above. In the case of even partial waves, this singularity 
is a pole, while for the odd partial waves the local frame 
transformation becomes singular by means of an infinite 
slope at e = iV. 

Next, let us consider the trace over the squared local 
frame transformations it/;/ [2^. Here, we observe the 
following behavior for approaching the threshold from 
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above 


lim ili/'(e) = ^ + lim 0(Ae) , with 

e\N Ae\,0 


l+l' 

cr’-E- 








— C(-^) 

p+1 ^ \ 2 


(38) 

(39) 


where the coefficients ^ are defined in Eq. ([S]) . In 
Eq. (l38l) we observe, that for the limit from above the 
elements of the matrix it smoothly approach the channel 
threshold irrespective of the particle exchange symmetry. 
However, the actual value attained does of course depend 
on the partial wave £, but there is no specific distinction 
between even and odd partial waves. On the contrary 
the limit from below exhibits a more intricate behavior 
and is given by 


lim it«'(e) = ^ + lim 

t/'N Ae\,0 


(S 


p—Q 




(iV-i)*^ 


(Ae)^ 


(40) 

where we observe that except from the constant value 
only terms contribute which have a fractional exponent 
in the energy dependence, implying non-analytic behav¬ 
ior when approaching the channel thresholds from below. 
We note that by comparing the local frame transforma¬ 
tion Uf,n and the elements of the closed channel coupling 
matrix il, their limiting behavior is just reversed, i.e. one 
approaches the limit in a singular manner from above, 
while regular from below and vice versa. 

Eor the case of even partial waves, it can be shown 
by exploiting general results on the Wigner 3j-symbols 
that the coefficients Cq^’^ ^ are always non-vanishing. This 
observation leads to the following asymptotic form 


lim ii«-(e) = 



lim (Ae) ^ 
Ae\,0 


(41) 


which means that the “du' ’s diverge at threshold as -^7=. 

Similar to the case of even partial waves we again ex¬ 
ploit general properties of the Wigner 3j-symbols to show 
for the odd values of i and £' that generally Cq^’^ ^ = 0, 

but c\’ ’ 0. This leads to the general from for the 

closed channel coupling, via 

lim il«'(e) = ^ -I- — f- lim (Ae)5 (42) 

Erom Eqs. dMli-dlll) we observe a fundamental differ¬ 
ence between even and odd partial waves. While the 
elements iiui for even partial waves are discontinuous 
across a channel threshold, the odd counterpart is con¬ 
tinuous. However, we note that the limit from above in 
the odd case approaches the value with an infinite 
slope, which is equivalent to say, that du' is not Lipshitz 
continuous across threshold. 



FIG. 11. (a) - (c) Transmission thresholds for the case ^ = 0. 
Close to threshold the transmission approaches unity with a 
universal slope independent of the scattering length. This 
universal behavior is present for all even partial waves in the 
^-SPWA. (d) - (f) Transmission around thresholds is shown 
for a variety of p-wave scattering lengths, demonstrating the 
dependence on the scattering length around threshold. 


The different threshold behavior for indistinguishable 
bosons and fermions, i.e. even and odd partial waves, is 
clearly observed in Eig. [TTJ where panels (a)-(c) depict 
the case for £ = 0 where, as expected from the preceding 
analysis, the transmission obtains a universal slope and 
value around the thresholds, regardless of the magnitude 
and the sign of the s-wave scattering length. This is due 
to the divergent behavior of il« and Ut,n across thresh¬ 
old. The trend exhibited in panels (a)-(c) is observed also 
for the first few even partial waves in £-SPWA. On the 
other hand, in panels (d)-(f) we illustrate the fermionic 
counterpart for £ = 1. Here, we observe a richer struc¬ 
ture of the threshold behavior due to the finite values the 
quantities iln and Uin acquire. This circumstance ren¬ 
ders the fermionic threshold behavior less universal in the 
sense that the transmission value obtained does depend 
on the scattering length. On the other side, this finite¬ 
ness observed for odd partial waves allows for a CIR even 
at threshold. Notably the same behavior is apparent in 
panels (d)-(f) for odd partial waves other than p-wave in 
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the £-SPWA picture. 


V. SUMMARY AND CONCLUSIONS 

The focus of our present study is on the scattering 
behavior of identical bosons and fermions as well as on 
distinguishable point-particles confined to a harmonic 
waveguide. The relevant coupling of different partial 
waves as well as the explicit energy dependence of the 
scattering properties are properly taken into account. 
The K matrix formalism established in [^, [13,113,113, 
is employed to obtain the relevant scattering observables. 
For higher partial wave interactions the explicit energy 
dependence of the scattering lengths is properly taken 
into account by including the free-space scattering re¬ 
sults of Gao |4^ . Using these results we are able 

to present fully analytical results including an adequate 
description for the interatomic scattering process. 

We consider the scattering process under the assump¬ 
tion of scale separation of the length associated with the 
interatomic and trapping potential, respectively. This as¬ 
sumption results in two regions of different symmetry, i.e. 
spherical close to the origin and cylindrical in the asymp¬ 
totic regime. The restriction that the total collision en¬ 
ergy lies below the threshold of the first excited transver¬ 
sal mode, as studied in previous works i HO, Hi was 
dropped and we thus allow for inelastic scattering de¬ 
scribing transversal (de-)excitation processes. The rela¬ 
tion between the quasi ID scattering amplitude and the 
physical K matrix found in was employed again to 
obtain analytical relations for the scattering observables. 
This formalism allows for a unified treatment of inelastic 
collisions within a harmonic waveguide for distinguish¬ 
able and indistinguishable particles. 

For the scattering of identical particles we investigated 
the transmission coefficient for up to four open channels 
where we find good agreement in the single partial waves 
approximation (SPWA) with the numerical results de¬ 
rived earlier [2l|. For the £-SPWA we also derived a quasi 
ID unitarity bound (Cl unitarity bound) which explains 
the influence of the open channels on the allowed trans¬ 
mission and transition coefficients around a CIR. The 
universal aspects of the Cl unitarity bound is demon¬ 
strated encapsulating all the relevant information of the 
interatomic interactions. However, the form of the Cl 
unitarity bound depends of the exchange particle symme¬ 
try. For the case of higher partial waves that are coupled 
we focus on the bosonic case where £ = 0 and 2 waves 
are coupled through the confinement. Here we studied 
the deviations from the d-SPWA to the case where the 
second partial wave is taken into account. In particu¬ 
lar we find that there is a region where the transmission 
when incident in the ground channel almost vanishes, 
while when incident in the first excited channel the par¬ 
ticles are non-interacting to the same degree. We studied 
this behavior for the free-space phase shifts as energy de¬ 
pendent quantities by using the analytic results of Gao 


[3^ on the scattering phase shift for potentials possessing 
a van der Waals tail, relating to recent experimental ob¬ 
servation ^^^Cs [4^, where d-wave shape-resonances were 
found. Furthermore, similar to the decoupling from the 
closed channels, described before by a vanishing element 
‘iXu we also find regions where the coupling between the 
partial waves vanishes, i.e. ilo2 vanishes. Over there, the 
d-SPWA is a good approximation to the transmission at 
a CIR. We note that also the corresponding fermionic 
case can be treated within the same formalism. For the 
case of distinguishable particles we derive resonance and 
transparency coefficients for the (dual) CIR. Here we ob¬ 
serve that the mechanism is way different from the case 
of indistinguishable particles where in particular the dual 
CIR was achieved by interference of partial waves coupled 
through the closed channels, while in the distinguishable 
case the coupling was accomplished by the open chan¬ 
nels. Finally we provide a brief description of the origin 
and the type of the appearing threshold singularities. 
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Appendix A: The physical K matrix 

In order to analytically invert the contribution of the 
closed channels to the physical K matrix, (1 — 
we employ the following method. First, we recognize, 
that the 3D K matrix, which is assumed to be diagonal, 
can be written as 

(Al) 

1=0 

here = tan 6i and denotes the i-th Cartesian basis 
vector. Using now the relation = C/^AT^^C/between 
the K matrices in the different regions in configuration 
space, where U denotes the local frame transformation 
[3514371 . we obtain 

^ AiU^eief U = ^ AAf, (A2) 

I 1^0 

where we have introduced the frame transformed basis 

(A3) 

For convenience we partition the frame transformed ba¬ 
sis into two parts f° and ff , corresponding to open and 
closed channels, respectively and also introducing the ab¬ 
breviations F° and to denote the dyads f° 0 f° and 
ff 0 fj, respectively. We note, that this definition gen¬ 
eralizes the open and closed K matrices attributed to a 


16 


specific partial wave, since the diagonal elements of F° 
and are given by 


A, F°u = 

A, F^u = Kc,i 


(A4) 

(AS) 


In particular we encounter the following frequently ap¬ 
pearing relation 

t^c r^lD \ _ 

\-^^oc,2-^ jk-^^co,l)'n,n' 

oo oo 

= E E 

r]—nc rj'—ric 

oo oo 

= AjA/ ^ ^ ^ ^ UinUirjUjrjUkr)^Ulr)'Uln' 

Tj—ric r}'—nc 

oo oo 

= A,AiU,nUln' UkrjUljj^ 

f]=nc •n=nc 

= A,AiikjiikliK)n,u', (A 6 ) 


(lAlOl) into the equation for the physical K matrix and 
using the relation ai/3ji = ctjPij one ends up with 

= detd -(k»°U) X (^'''« + 

- lAeAg, {iiti,F°u + iI«F?,^, - iI«-(F^,, + F^,^))), 

(A13) 

representing the physical K matrix interacting with two 
partial waves i. and . We note, that this formula only 
holds for indistinguishable particles, i.e. both partial 
waves have to be either even or odd. In addition, we 
readily observe that the K matrix given in Eq. (IA13I) is 
real and symmetric as expected. 


Appendix B: The resonance and transparency 
coefficients 


where n, n! range within the open channels and the ily’s 
are the coupling elements derived earlier [2^. These en¬ 
ergy dependent elements are defined according to 

OO 

ii«' = ^ UlnUVn. (A7) 

n—Tio 

which are given in Eqs. o, m and (| 6 ]) in closed form. 
In a similar fashion to Eq. (lAfil) . the following relation is 
derived 

= A,A, il.,F° (A 8 ) 

Appearing higher order products of the closed channel K 
matrices are readily shown to satisfy 


This section is considered as a brief summary of notions 
used here, which were introduced in in order to keep 
this presentation as self-contained as possible. 

We remind that the formation of a CIR is described 
by a diverging physical K matrix, i.e. the roots of 
det(l — iKcc)- As it can be seen clearly in Eq. (IA13I) this 
divergence can only be achieved by a vanishing denom¬ 
inator. Equivalently, this can be expressed by a diver¬ 
gence of the couplings, given in Eqs. (lAllI) and (IA12I) . 
Explicitly this means that a £-wave CIR occurs in the 
SPWA when ai diverges, and similar, a CIR oc¬ 

curs in the presence of Awave interactions, when di¬ 
verges. Parameterizing this divergence of in terms of 
the scattering length and energy, yields 


= A,A,Afc (A9) 

The relations (A 6 ) to (A 8 ) turn out to be very useful 
when actually carrying out the analytical inversion of the 
matrix (1 — Similar to the procedure in [1^, the 

inversion is done by first recognizing, that the expansion 
of the ID K matrix given in Eq. (IA2I) is written as a sum 
over dyads, i.e. rank one matrices and then repeatedly 
applying the Sherman-Morrison formula [1^ . The result 
of this procedure yields 

(1 - )-i = 1 + a,(l + 

+ (F)). + F^%), (AlO) 


ai{e) 

RC^(e) 


RC,(e) 

-1 


2 y^ e + 1/2 


iil«(e) ’ 


(Bl) 

(B2) 


and similar when the coupling between two partial waves 
has to be taken into account 


= RCf/^(e) 

RC.^Ke) = 


2 -^/e -|- 1/2 


X 




(B3) 

(B4) 


where the coefficients a and /? are given by 
iAi 

Oti = - 

1 - iAiiiu 

R-. = __ 

l-*A,(%-a.Uf^.) 


(All) 

(A12) 


These coefficients play a similar role as the couplings gi 
and 52 defined in |24| . Inserting now the result of Eq. 


The so-called dual CIR, where the total transmission 
becomes unity, is obtained within the same framework 
by a vanishing numerator of the physical K matrix and 
is due to the matrix nature of K in the case of multiple 
open channels only expressible in terms of transparency 
coefficients Tonly in the case of a single open channel. 
However, by analogous arguments as for the resonance 
coefficients RC, introduced above, the parametrization of 
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a dual CIR is then given by 


ai'ie) 

TQ.,,(e) 


TQ.,(e) 

1 

2 -^/ e + 1/2 


where we note that due to the needed destructive inter¬ 
ference, the dual CIR is only possible when more than 
(B5) one partial waves are taken into account. 


21 ’+ 1 
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+ iiuUlo - 2ilu,UioUe>o) '' 
(B6) 
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